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Inelastic  dipole  scattering  of  electrons  from  surface  vibrational 
resonances  of  clean  Cu(110)  and  Ni(110)  surfaces  has  been  observed. 
The  resonances  are  shown  by  surface  lattice  dynamics  calculations 
to  be  a  characteristic  signature  of  the  (110)  surface  of  fee  metals. 
The  calculated  spectral  energies  and  shapes  are  in  excellent  agree¬ 
ment  with  experiments  and  shows  that  there  are  no  significant 
changes  in  the  surface  force  constants  from  the  bulk  values. 
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Vibrational  spectroscopy  of  surface  vibrational  modes  on  clean  and  adsor¬ 
bate  covered  metal  surfaces  is  a  viable  method  for  obtaining  valuable  knowledge 
about  the  structure  and  the  interatomic  forces  at  the  surface.  1,2  For  a  few  close 
packed  surfaces  it  has  recently  been  possible  by  electron-energy-loss  spectroscopy 
(EELS)  3  and  inelastic  He  scattering  4  to  observe  surface  phonon  dispersion  rela¬ 
tions  ,  modes  attributed  to  surface  resonances  5  and  dipole  excited  bulk  modes.  6 
In  this  paper  we  present  observations  by  EELS  of  dipole  active  surface  resonances 
on  the  (110)  surfaces  of  Cu  and  Ni.  Calculations  using  surface  lattice  dynamics 
are  in  excellent  agreement  with  measured  spectra  and  show  that  these  resonances 
have  their  physical  origin  in  the  geometrical  structure  of  the  (110)  surfaces  of 
face  centered  cubic  crystals.  In  particular  it  is  found  that  there  are  no  significant 
changes  in  the  surface  force  constants  from  the  bulk  values. 

The  EEL  spectra  were  recorded  in  a  multitechnique  ultrahigh  vacuum  sys¬ 
tem  with  a  base  pressure  of  4xl0~n  Torr.  7  The  EEL  spectrometer  consists  of  a 
monochromator  and  analyzer  of  the  double  pass  127  ’  cylindrical  electrostatic 
deflector  type.  The  spectra  were  recorded  in  the  specular  direction  with  an 
incident  angle  of  60  ’  with  respect  to  the  crystal  normal.  The  scattering  plane 
containing  the  incident  and  scattered  electrons  is  defined  by  the  surface  normal 
and  the  (110)  crystal  direction.  The  angular  acceptance  of  the  analyzer  is  1.8  *  at 
full-width-half-maximum  (FWHM).  Impact  energies  of  3.2  and  4.3  eV  were  used 
for  Cu  and  Ni,  respectively.  Sample  cleaning  and  characterization  were  done  by 
standard  procedures.  7 

The  recorded  spectra  in  Fig.l  display  sharp  spectral  features  which  are  only 
slightly  broader  than  the  clastic  peak.  The  observed  losses  show  the  expected 
reduction  in  intensity  from  the  Bose-Einstein  factor  on  the  gain  side.  Off  specular 
measurements  show  that  the  inelastic  scattering  cross  section  follows  closely  those 
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of  dipole  active  modes.  The  loss  peaks  occur  at  energies  20  and  24  meV,  which 
scale  with  and  lie  well  below  the  maximum  bulk  phonon  energies  at  30  8  and  37 
meV  9  for  Cu  and  Ni,  respectively.  This  fact  suggests  that  these  vibrational 
features  are  derived  from  the  bulk  phonons  of  the  crystals.  The  electric  field 
from  the  scattering  electron  can  also  excite  electron  hole  pairs.  10  In  general  this 
mechanism  gives  rise  to  an  additive  contribution  to  the  spectrum  which  decreases 
monotonically  with  frequency  and  accordingly  can  not  be  responsible  for  the 
observed  sharp  loss  features. 

The  inelastic  dipole  scattering  from  phonons  at  metal  surfaces  is  readily 
related  to  surface  lattice  dynamics  of  dipole  active  displacements.  The  lattice 
dynamics  can  be  done  analytically  by  a  transfer  matrix  technique  for  a  semi¬ 
infinite  surface.  11 

The  probability  distribution  P[lt,k  )  for  an  electron  with  incident  wavevec- 
tor  k  to  be  detected  in  the  direction  it  with  an  energy  loss  Hu  and  for  small 
momentum  transfers  ft7||  —  ^\\  ~  ^||  parallel  to  the  surface  is  given  by  1,12 


8  *' 
ftjr(efl0)2cosar  k 


n 

(9  |f+7..2)3 


(1) 


where  a  is  the  angle  of  incidence,  a0  the  Bohr  radius,  z  the  surface  normal  and 
qt  =  |  kt  +kt'  | .  The  prefactor,  i.e.  the  dipole  matrix  element,  to  the  spectral 
function  S(7||,^)  in  Eq.(l)  is  sharply  peaked  close  to  7||  =  0.  In  this  limit  the 
spectral  function  is  given  by  a  temporal  Fourier  transform  of  a  correlation  func¬ 
tion  for  the  normal  component  of  the  dipole  moment  /*  of  the  charge  fluctuations 
in  the  metal. 

For  the  high  symmetry  surfaces  only  a  rigid  displacement  ti(  of  a  layer  l 
normal  to  the  surface  can  give  rise  to  a  dipole  moment  ft.  This  means  that  the 
incident  electron  will  essentially  excite  only  longitudinal  phonons  propagating 


along  the  direction  normal  to  the  surface  by  dipole  scattering  for  the  (110),  (111) 
and  (100)  surfaces  of  Cu  and  Ni. 


The  effective  charge  relates  the  dipole  moment  /i  to  the  displacements  u, 

of  the  substrate  layers  through  /j=Ve;*U{.  Recently  it  has  been  shown  that  the 

l 

incomplete  screening  by  the  conduction  electrons  in  the  surface  region  can  give 
rise  to  sufficiently  large  et*  so  that  the  dipole  coupling  to  the  bulk  phonons  of  the 
(100)  surfaces  of  Cu  and  Ni  is  observable.  8  Here  the  same  model  is  used  for  the 
effective  charges  where  t[  = and  e/*  =  0  for  />3.  Note  that 
the  sum  of  the  effective  charges  has  to  be  zero  since  a  rigid  translation  of  the 
whole  crystal  normal  to  the  surface  does  not  change  the  dipole  moment. 

In  the  harmonic  approximation  the  spectral  function  5(y||=0,u;)  can  be 
directly  related  to  a  phonon  density  of  states  jf(w,{n,'})  projected  onto  a  normal¬ 
ized  effective  charge  field  nt'  =  t\f  et"t  through  12 


5,  Tl-M  =  (».•)) 


(2) 


Here  M  is  the  mass  of  a  substrate  atom,  n  (w)  the  Bose-  Einstein  distribution  fac¬ 
tor,  a  the  area  of  the  surface  primitive  cell  and  et*t  the  modulus  of  the  effective 
charge  field  {elr}  ,  i.e.  in  the  adopted  model  e(*t  =  v/2e  *  ,  and 


=  E  I  E ni'uKi  I  > 

K  t 


(3) 


where  uk/  is  the  rigid  displacement  of  a  substrate  layer  /  for  an  eigenmode  k 
with  energy  of  the  dynamical  matrix  D  . 

A  prerequisite  for  a  meaningful  calculation  of  jj(u,’,{n/})  is  a  good  descrip¬ 
tion  of  the  force  constants  with  respect  to  the  rigid  displacements  of  the  sub¬ 
strate  layers  for  Cu  and  Ni.  As  shown  in  Fig.  2,  a  simple  oue-parameter  nearest 


-  5- 


neighbor  central  force  constant  model  between  the  atoms  gives  a  good  account  of 
the  measured  dispersion  of  phonons  propagating  along  the  [110]  direction  for  both 
Cu  and  Ni.  8,9  This  dispersion  relation  shows  a  local  maximum  well  within  the 
Brillouin  zone  boundary  in  contrast  to  the  calculated  dispersion  for  phonons 
along  the  other  two  directions.  This  behavior  is  a  direct  consequence  of  the  fact 
that  an  atom  in  a  metal  layer  normal  to  the  [1 10]  direction  has  nearest  neighbor¬ 
ing  atoms  in  both  the  first  and  second  neighboring  layers  while  it  has  only 
nearest  neighbors  in  the  first  neighboring  layers  for  layers  normal  to  the  [ill]  or 
the  [100]  directions.  Thus  the  eigenvalue  problem  for  the  dynamics  of  the  pho¬ 
nons  in  the  [110]  direction  will  have  up  to  second  nearest  interlayer  force  con¬ 
stants  ; 


*«l  -  «<- 2  -  1  -  «/+l  -  «<+  2) 


(4) 


Here  the  force  constant  k  is  directly  related  to  the  frequency  um!U  by 
k  (M  =  Wm„/8  .  Note  that  the  second  nearest  interlayer  force  constant  is  as 
strong  as  the  first  interlayer  force  constant.  This  causes  the  restoring  force  for  a 
displacement  field  to  be  larger  at  the  reduced  wave  vector  f=0.5  than  at  the  zone 
boundary  f*=1.0.  For  the  other  two  directions  the  geometrical  configuration  of 
the  atoms  causes  the  interlayer  force  constants  to  be  twice  as  large  and  to  couple 
only  to  the  nearest  neighboring  layers.  This  causes  the  restoring  force  to  be 
monotonic  with  f  and  to  be  larger  than  the  corresponding  forces  for  the  [110] 
direction  at  j=l.0. 

The  loss  of  coordination  for  the  atoms  at  the  surface  is  the  only  effect  taken 
into  consideration  in  modelling  the  inter-  and  intralaycr  surface  force  constants. 
In  the  case  of  the  (110)  surface  this  gives  rise  to  the  following  modification  of  Eq. 
(•I)  at  the  two  outermost  layers 


All  the  stationary  points  in  the  dispersion  relations  for  the  phonons  in  Fig.  2 


give  rise  to  divergent  van  Hove  singularities  in  the  density  of  states  projected  on 


the  bulk  layers  as  shown  in  Fig.  3  .  These  divergences  disappear  in  the  projec¬ 


tion  on  the  surface  layers,  due  to  the  destructive  interference  between  the 


incident  and  reflected  phonons  from  the  surface  at  these  energies.  For  the  (110) 
surface,  there  is  a  pseudo  band  gap  (c.f.  Fig.3  upper  panel)  for  0  <  u  <  --’rnax/^ 


,  and  for  every  energy  in  this  range  there  exists  an  evanescent  wave  localized  to 


the  surface.  At  a  certain  energy  depending  on  the  surface  force  constants  this 


wave  will  have  minimal  coupling  through  Eq.(5)  to  the  bulk  phonons  and  will 
show  up  as  a  resonance  as  seen  in  Fig.  3.  For  the  (111)  and  (100)  surfaces  there 
are  no  distinct  spectral  features  in  the  corresponding  density  of  states. 


A  detailed  comparison  of  the  calculated  resonance  with  the  experimental 
spectra  has  been  done  by  including  the  frequency  dependent  prefactors  appearing 
in  Eqs.  (1)  and  (4),  and  by  integration  of  the  inelastic  scattering  probability  over 
the  acceptance  aperture  of  the  detector.  The  resulting  electron  energy  loss  func¬ 
tion,  I(u>)  (inset  in  Fig.  1),  shows  a  rather  sharp  loss  feature  which  agrees  closely 
with  the  measured  loss  peak  and  is  only  2  meV  broader  than  the  instrumental 
resolution.  This  agreement,  equally  good  for  the  corresponding  Cu  surface, 
implies  that  the  changes  in  the  surface  interlayer  force  constants  from  the  bulk 
values  are  no  more  than  ±15?o  .  The  measured  spectral  loss  intensities  in  Fig.  1 
are  reproduced  with  the  effective  charges  t  *  =  0.024±0.003  and  0.028±0.003  for 


Cu  and  Ni,  respectively,  which  are  only  slightly  larger  than  the  value  0.016 


obtained  for  Cu  atoms  in  the  surface  region  ®. 


/v'.-V „--w.  ..  v  .• 


In  summary  surface  resonances  of  longitudinal  bulk  phonons  have  been 
observed  on  the  (110)  surfaces  of  Cu  and  Ni.  These  sharp  resonances  are  shown 
to  be  a  general  feature  of  the  (110)  surface  of  face  centered  crystals  with  peak 
energies  scaling  with  the  maximum  bulk  phonon  energies.  Detailed  comparison  of 
measured  spectra  with  results  from  surface  lattice  dynamics  shows  that  there  are 
no  significant  stiffening  or  softening  of  the  surface  force  constants. 
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Council  is  gratefully  acknowledged. 


References 


t  On  leave  from:  Institute  of  Theoretical  Physics,  Chalmers  University  of  Tech 
nology,  S-412  96  Goteborg,  Sweden. 

Present  address:  Materials  and  Molecular  Research  Division,  Lawrence  Berkeley 
Laboratory,  Berkeley,  CA  94720. 

1.  H.  Ibach  and  D.L.  Mills,  Electron-Energy-Loss  Spectroscopy  and  Surface 
Vibrations  (Academic,  New  York,  1982). 

2.  See  articles  in  Proceedings  of  the  Third  International  Conference  on  Vibra¬ 
tions  at  Surfaces,  Asilomar,  California,  edited  by  C.R.  Bruudle  and  H. 
Morawitz,  J.  Electron  Spectrosc.  Relat.  Phenom.  29  and  30  ,  (1983). 

3.  S.  Lehw'ald,  J.M.  Szeftel,  H.  Ibach,  T.S.  Rahman,  and  D.L.  Mills,  Phys.  Rev. 
Lett.  60 , 518  (19S3) 

4.  R.B.  Doak,  U.  Harten  and  J.P.  Toennies,  Phys.  Rev.  Lett.  61  ,  578  (1983); 
M.  Cates  and  D.R.  Miller,  J.  Electron  Spectrosc.  Relat.  Phenom.  30  ,  157 
(1983) 

5.  V.  Bortolani,  A.  Franchini,  F.  Nizzoli,  and  G.  Santoro,  Phys.  Rev.  Lett.  62 
,  429  (1984). 

6.  S.  Andersson,  B.N.J.  Persson,  M.  Persson,  and  N.D.  Lang,  Phys.  Rev.  Lett. 
62 , 2073  (1984). 

7.  J.A.  Stroscio,  S.R.  Bare,  and  W.  Ho,  Surf.  Sci.,  in  press. 

8.  E.C.  Svensson,  B.N.  Brockhouse,  and  J.M.  Rowe,  Phys.  Rev.  166  ,619 
(1967). 

9.  R.J.  Birgcnau,  J.  Cordes,  G.  Dolling,  and  A.D.B  Woods,  Phys.  Rev.  136  , 


A 1359  (1964). 


Figure  Captions 

1.  Electron  energy  loss  spectra  of  clean  Cu(110)  (upper  panel)  and  Ni(UO) 
(lower  panel)  surfaces.  The  spectra  were  recorded  at  300  K  in  the  specular 
direction.  Sharp  spectral  features  are  observed  both  on  the  loss  and  the  gain 
side  at  20  and  21  meV  for  Cu  and  Ni,  respectively  .The  inset  shows  the  cal¬ 
culated  electron  energy  loss  function  I(u)  for  the  clean  Ni(110)  surface  nor¬ 
malized  to  the  maximum  elastic  intensity.  The  instrumental  resolution  is 
accounted  for  by  Gaussian  broadening  with  4  meV  at  full-width-half¬ 
maximum  (FWHM).  The  effective  charges  t  *  and  -t  ‘  of  the  two  outer¬ 
most  surface  layers  have  been  adjusted  to  reproduce  the  measured  spectral 
loss  intensities. 

2.  Bulk  phonon  dispersion  relations  in  reduced  units  along  different  major  cry¬ 
stal  directions  for  Cu  and  Ni.  The  calculated  dispersion  for  longitudinal  pho¬ 
nons  propagating  along  the  [110]  direction  (solid  curve)  are  compared  with 
measurements  by  inelastic  neutron  scattering  8,0  The  experimental  results 
have  been  scaled  with  the  frequency  u;mM  where  wmaix/v/2  is  the  measured 
frequency  at  the  Brillouin  zone  boundary,  urmix  =  29.7  and  36.7  meV  for  Cu 
and  Ni,  respectively.  The  dashed  curve  is  the  calculated  dispersion  for  the 
longitudinal  bulk  phonon  along  the  crystal  directions  (100)  and  (ill).  The 
distance  between  the  crystal  layers  is  d. 

3.  Projected  density  of  states  for  the  longitudinal  phonons  in  Fig.  2.  They 
have  been  calculated  for  a  projection  on  the  relative  displacement  of  two 
neighboring  bulk  layers  (upper  panel)  and  the  two  outermost  surface  layers 
(lower  panel).  The  resonance  in  the  density  of  states  projected  on  the  surface 
layers  for  the  phonons  propagating  along  the  (110)  direction  (solid  curve  in 
the  lower  panel)  peaks  at  19.4  and  21.0  meV  for  Cu  and  Ni,  respectively. 
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